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ABSTRACT 


This  thesis  evaluates  the  protection  provided  by  towed-decoys  deployed  by  an 
aircraft  during  an  engagement  against  an  anti-air  missile  equipped  with  a  monopulse 
seeker.  The  research  emphasizes  the  use  of  passive  decoys.  Many  of  the  operational 
parameters  required  before  the  deployment  of  a  towed-decoy  are  investigated,  including 
the  strength  of  reflection,  the  tether  length,  and  the  direction  of  release,  under  different 
missile  incoming  directions. 

The  thesis  evaluates  two  reflection  cases.  One  is  for  both  the  target  and  the  decoy 
being  point  reflectors.  The  other  is  for  an  extended  target  and  a  point  reflected  decoy.  The 
protection  envelopes  for  different  engagement  scenarios  of  each  reflection  case  is 


obtained. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

n.  MONOPULSE  ANGLE  ESTIMATION .  3 

in.  SIMULATION  ARCHITECTURES . 7 

IV.  DECOY  EFFECTS  FOR  A  POINT  REFLECTION  TARGET  9 

A.  LAYOUT  GEOMETRY . 9 

B.  EFFECTS  OF  REFLECTION  RATIO .  10 

C.  CHANGE  OF  MISS  DISTANCES  DUE  TO  DIFFERENT 

MISSILE  ATTACK  ANGLES . H 

D.  CHANGE  OF  MISS  DISTANCES  DUE  TO  DIFFERENT 

DECOY  DEPLOYMENT  ANGLES . 13 

V  DECOY  EFFECTS  FOR  AN  EXTENDED  TARGET . 15 

A.  EFFECTS  OF  REFLECTION  RATIO . 16 

B.  CHANGE  OF  MISS  DISTANCE  DUE  TO  DIFFERENT 

MISSILE  ATTACK  ANGLES . 16 

C.  CHANGE  OF  MISS  DISTANCE  DUE  TO  DIFFERENT 

DECOY  DEPLOYMENT  ANGLES . 17 

VI.  CONCLUSIONS .  19 

VII.  FUTURE  STUDIES . 21 

APPENDIX  A  -  THESIS  MAIN  PROGRAM . 23 


vii 


APPENDIX  B  -  SEEKER  MODEL 


29 


LIST  OF  REFERENCES . 31 

BIBLIOGRAPHY . 33 

INITIAL  DISTRIBUTION . 35 


LIST  OF  FIGURES 


1.  Scenario  Geometry . ^ 

2.  Radar  Tracking  Axes . ^ 

3.  Antenna  Patterns  and  Error  Voltage . 4 

4.  Angle  Tracking  Block  Diagram . 5 

5.  Simulation  Block  Diagram . ^ 

6.  Layout  Geometry . ^ 

7.  Effect  of  Reflection  Ratio  on  Miss  Distance, 

Decoy  at  150®  and  Missile  from  180® . 

8.  Effect  of  Reflection  Ratio  on  Miss  Distance, 

Decoy  at  180®  and  Missile  from  90® . 11 

9.  Protection  Envelope  of  Different  Missile  Attack  Angle 

for  <Jo/Ox  =  . . 

10.  Protection  Envelope  of  Different  Missile  Attack  Angles 

for  Oj/Ot  “  . . 

11.  Protection  Envelope  of  Different  Decoy  Deployment  Angles 

for  OyJCj  =  1.5 . 14 

12.  Extended  Target  Reflection . 15 

13.  Effect  of  Reflection  Ratio  on  Miss  Distance 

for  Extended  Target . 16 

14.  Protection  Envelope  of  Different  Missile  Attack  Angles 

for  Extended  Target . 17 


IX 


15.  Protection  Envelope  of  Different  Decoy  Deployment  Angles 
for  Extended  Target . 


18 
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I.  mXRODUCTION 


Monopulse,  or  simultaneous-lobing  tracking,  obtains  an  angle  error  estimate  from 
a  single  pulse  return.  Unlike  sequential-lobing  systems,  monopulse  radar  does  not  require 
a  long  integration  time  or  demodulation  process,  resulting  in  an  immunity  to  noise 
jamming  effects.  This  inherent  ECCM  against  many  on-board  jamming  techniques  explains 
why  so  many  missiles  have  adopted  monopulse  trackers  for  their  guidance  seekers.  Also 
many  sequential-lobing  trackers  have  been  modified  using  monopulse  techniques  to  reduce 
their  vulnerability  to  jamming. 

From  a  self-protection  point  of  view,  a  fighter  pilot  needs  an  effective  ECM 
method  to  counter  missiles  with  monopulse  seekers.  Since  on-board  equipment  has  limited 
effectiveness,  expendables,  especially  decoys,  become  a  reasonable  choice.  However, 
modem  missiles  can  distinguish  the  difference  in  Doppler  frequencies  between  aircraft  and 
conventionally  free-falling  decoys  using  MTI  techniques  leading  to  a  preference  for 
towed-decoys. 

An  initial  evaluation  leads  to  the  conclusion  that  a  towed-decoy  provides  a  counter 
to  a  monopulse  seeker.  However,  its  effectiveness  is  always  in  question.  A  critical 
question  concerns  the  protection  it  can  provide  under  different  missile  attack  angles. 
Further,  there  are  some  parameters  that  must  be  decided  prior  to  deployment  of  a  decoy. 
For  example,  the  length  of  the  tether  line,  the  strength  of  the  reflection  from  the  decoy, 
and  the  direction  of  release,  to  name  just  a  few.  This  thesis  investigates  all  the  above 
questions  using  a  computer-based  simulation  based  upon  the  scenario  shown  in  Figure  1. 
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Figure  1.  Scenario  Geometry 


Description  of  Parameters: 

X„  Yj  =  Inertial  reference  axes 

Rj^„  Rp,  Rjj  =  Missile,  target,  and  decoy  position  vectors  with  respect  to  inertial 

frame 

Vg  =  Missile  velocity  vector 

=  Position  vector  of  missile  to  target 
Rnpd  =  Position  vector  of  missile  to  decoy 
L  =  Tether  line  length 
y  -  Decoy  deployment  angle 
a  =  Subtended  angle  between  target  and  decoy 

67,  ©M,  6d  “  Orientations  of  target,  missile,  and  decoy  with  respect  to  inertial  frame 
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n.  MONOPULSE  ANGLE  ESTIMATION 


There  are  several  types  of  monopulse  angle  sensors  that  are  categorized  by  their 
sensing  and  detection  methods.  The  most  commonly  used  and  preferred  one  is  the 
amplitude  sensing,  sum-difference  detection  type.  This  is  the  kind  of  sensor  investigated  in 
this  thesis. 

The  general  radar  seeker  model  is  shown  in  Figure  2. 


In  Figure  2,  the  reference  axis  is  the  missile  body  roll  axis  while  0^  is  the  angular 
tracking  error  defined  as  the  angle  between  the  seeker  tracking  axis  and  the  line-of-sight. 
Assuming  Gaussian  distributed  antenna  patterns: 

G(0)  =  Goexp{-k(^)2}  (1) 

where,  03  is  the  antenna  beamwidth,  0^  is  the  squint  angle,  and  k  =  2.776,  results  in 
received  voltages  given  by: 
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Amplitude 


Vi  =  SGoexp{-k(^)^)  and  V2  =  SGoexp{-k(^)^}  (2) 

and  S  is  the  strength  of  the  target  skin  return.  The  receiving  patterns  are  shown  in  Figure 
3-a. 


Figure  3-a.  Gaussian  Antenna  Patterns 

3-b.  Sum  Beam  and  Difference  Beam  Patterns 
3-c.  Error  Voltage 


If  the  IF  gains  of  both  sum  and  difference  channels  are  matched  (i.e.,  Kj.  =  K^), 
then,  the  outputs  of  the  sum  and  difference  channels  are: 

L  =  V1+V2  and  A  =  vi -V2  (3) 


Now,  the  output  of  the  phase  detector  is  the  error  voltage,  e, 
e  =  ReaKViA^i}  (4) 
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The  error  voltage,  or  monopulse  ratio,  is  proportional  to  the  angular  tracking 
error,  0^.  It  serves  as  an  input  to  the  seeker  angle  servo  which  acts  to  lock  the  tracking 
axis  onto  the  target  by  nullifying  the  tracking  error.  The  sum  and  difference  receiving 
channels  signals  and  resulting  error  voltage  are  shown  in  Figure  3-b  and  3-c.  The  whole 
process  is  illustrated  in  the  block  diagram  of  Figure  4. 


Figure  4.  Angle  Tracking  Block  Diagram 


As  shown  in  Eq.  4,  the  error  voltage  is  the  real  part  of  the  normalized  (with 
respect  to  the  sum  signal)  difference  signal.  No  matter  how  strong  an  on-board  jamming 
signal  is,  it  is  always  spatially  coincident  with  the  target  signal.  Thus,  the  error  voltage 
does  not  change  with  on-board  jamming.  That  is  why  an  off-board  source  is  needed  to 
perform  angular  deception  against  a  monopulse  seeker. 

Double-Source  Problem 

For  a  double-source  geometry  shown  as  Figure  1  of  Chapter  I,  the  received 
voltage  from  each  beam  is: 

vi  =  vt,  +vdi  =  SGoexp{-k(^^^)^}  +  JGoexp  |-k(-^^)^|ej'''  (5) 

V2  =  VT,  +VD,  =  SGcexp{-k(2^)"l  +  JGoexp{-k(2^)")ei''  (6) 
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where  \|/  is  the  phase  relation  between  target  and  decoy  returns.  This  is  given  by: 

V  =  (7) 

where  5L  is  the  path  difference.  Then,  the  sum  and  difference  channels  have  voltages: 

Z  =  vi  +  V2  and  A  =  vi  -  va  (8) 

Again,  the  indicated  angle  can  be  derived  from  the  error  voltage: 

e  =  Real{A/Z}  (9) 

The  introduction  of  the  second  source  results  in  distortion  of  the  real  part  of  the 
monopulse  ratio.  The  tracking  axis  does  not  point  at  the  target  any  longer.  Instead,  it 
generally  tracks  the  power  centroid  of  the  two  sources  until  they  can  be  resolved.  When 
the  seeker  starts  to  resolve  the  two  sources,  it  tracks  the  stronger  source.  The  resolving 
condition  occurs  when  the  subtended  angle  is  greater  than  the  seeker's  beamwidth;  this 
distance  is; 


s  <(j)/tan(f)  =  L/(2tan^) 


(10) 
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m.  SIMULATION  ARCHITECTURES 


The  simulation  is  explained  using  the  block  diagram  shown  in  Figure  5.  Using 
missile/target  geometry,  the  aspect  angles  of  the  target  and  decoy  are  provided  as  inputs 
to  the  seeker  model.  This  allows  the  tracking  error  to  be  estimated.  The  seeker's  tracking 
axis  is  pointed  at  the  line-of-sight  between  the  missile  and  the  composite  target/decoy 
return.  The  lateral  acceleration  required  to  turn  the  missile  is  simply 

acMD  =N»6a*Vb 

where  N  is  the  navigation  constant.  This  steering  command  goes  into  the  actuator  model 
to  drive  the  deflection  of  the  missile  control  surface  and  change  the  missile  attitude 
towards  the  target. 

To  simplify  the  analysis  several  assumptions  are  made.  First,  the  simulation  is  limited 
to  an  azimuth-plane  only  rather  than  to  a  three-dimensional  geometry.  This  simplifies  the 
simulation  and  allows  results  to  be  obtained.  Second,  the  missile  dynamics  are  assumed  to 
have  a  1  Hz  response  while  the  actuator  is  limited  to  lOG  acceleration.  Third,  the 

monopulse  seeker  uses  amplitude  sensing  with  sum-difference  detection.  Further  it  is 
assumed  that  the  seeker  gimbal  has  20  Hz  response  with  angle  limits  of  ±45°.  This  high 
dynamic  seeker  gimbal  performance  avoids  dynamic  coupling  with  the  missile.  The 

antenna  beamwidth  (63)  is  approximately  1 1.5®. 

The  measure  of  effectiveness  used  is  the  miss  distance  (i.e.,  the  range  between 
missile  and  target).  It  is  set  such  that  the  missile  will  detonate  if  it  is  within  5  meters  range 
to  the  target  (either  the  aircraft  or  the  decoy).  There  are  times  when  the  missile  does  not 
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impact  on  either  of  the  targets.  Under  these  circumstances,  the  miss  distance  is  recorded 
as  the  shortest  missile-to-aircraft  range  (  i.e.,  minimum  during  the  whole  simulation). 
Therefore,  these  recorded  miss  distances  become  the  minimum  protection  envelope  that 
the  decoy  can  provide. 


Figure  5.  Simulation  Block  Diagram 
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IV.  DECOY  EFFECTS  FOR  A  POINT  REFLECTION  TARGET 


A.  LAYOUT  GEOMETRY 

Before  discussing  the  effects  of  different  decoy  deployment  conditions,  the 
reference  coordinate  system  is  first  defined.  There  are  two  angles  of  interest  in  this  thesis. 
One  is  the  decoy  deployment  angle  while  the  other  is  the  initial  missile  attack  angle.  Let 
the  angle  P  be  defined  as  shown  in  Figure  6.  In  Figure  6-a,  P  =  y,  is  used  to  analyze  decoy 
deployment  angle  effects;  whereas  in  Figure  6-b,  P  =  0^,  is  used  to  analyze  initial  missile 
attack  direction  effects.  The  missile/target  engagement  situation  can  be  categorized  as: 

(i) .  For  P  =  45”  ~  135"  or  -45"  ~  -135",  the  target  aircraft  undergoes  a  side-attack 
situation. 

(ii) .  For  P  =  135"  ~  225",  the  target  aircraft  undergoes  a  head-on  attack  situation. 

(iii) .  For  other  P,  the  target  aircraft  undergoes  a  tail-pursuit  situation. 
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B.  EFFECTS  OF  REFLECTION  RATIO 


Figure  7  shows  the  effects  on  miss  distance  resulting  from  three  different  tether 
lengths  as  a  fiinction  of  the  relative  decoy  to  target  ratios.  It  shows  that  the  longer  the 
tether  line  is,  the  higher  the  reflection  ratio  (aj/aT.)  required  to  effectively  capture  the 
seeker's  angle  tracking  loop.  Moreover,  when  the  reflection  ratio  reached  1.5,  the  effects 
on  miss  distance  stabilized.  At  that  point,  the  missile  will  home  onto  the  decoy. 


Figure  7.  Effect  of  Reflection  Ratio  on  Miss  Distance, 
Decoy  at  ISO”  and  Missile  from  180' 


Figure  8  shows  a  side-attack  situation.  The  decoy  is  more  effective  with  this 
operation,  and  =  1.25  is  sufficient  to  protect  the  aircraft. 


10 


O 

c 

OS 

lUU 

oS 

b 

80 

u> 

</) 

2 

60 

40 

20 

/  L=1C 

)0m 

/  /  L=75 

/  / 

/  / 

/  ✓  _ 

m 

/ 

/  / 

'  /  L=50 

m’"’ 

— 

17 

/li 

/  P 

y  / 

y  f 

V 

Figure  8.  Effect  of  Reflection  Ratio  on  Miss  Distance, 
Decoy  at  180®  and  Missile  from  90® 


C.  CHANGE  OF  MISS  DISTANCES  DUE  TO  DIFFERENT  MISSILE 
ATTACK  ANGLES 

Figure  9  shows  that  for  a  reflection  ratio  of  1  with  the  decoy  deployed  at  1 80®,  the 
protection  from  tail-approaching  missiles  stays  almost  constant  independent  of  tether 
lengths.  The  protection  from  side-attack  missile  reduces  as  the  tether  length  increases. 
Note  also  that  when  the  tether  length  increases,  null  regions  become  more  obvious  and  the 
protection  regions  get  smaller.  This,  again,  shows  that  for  longer  tether  lengths,  a  higher 
reflection  ratio  is  required  to  capture  the  angle-tracking  loop  of  the  seeker. 
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This  point  is  illustrated  by  Figure  10  where  for  a  reflection  ratio  of  1.5  a  wider 
protection  envelope  is  provided.  Effective  protection  is  provided  when  encountering  a 
tail-attack  missile,  with  the  decoy  parameters  set  to  =1.5  and  L  =  50m;  whereas 
when  a  side-attack  missile  is  encountered,  =1.5  and  L  =  75m  is  recommended.  The 
reason  that  L  =  100m  isn't  preferred  is  due  to  the  many  nulls  in  the  protection  pattern, 
which  means  that  in  some  cases  the  missile  may  lock-on  to  the  aircraft.  Also,  long  tether 
lengths  are  not  convenient  if  evasive  maneuvers  are  required. 
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Figure  10.  Protection  Envelope  of  Different  Missile  Attack  Angles 
for  CyjGj  =  1.5 


D.  CHANGE  OF  MISS  DISTANCES  DUE  TO  DIFFERENT  DECOY 
DEPLOYMENT  ANGLES 

The  following  discusses  the  parameters  associated  with  decoy  deployment  for 
cJCj  =1.5  with  the  missile  attacking  from  180".  Figure  11  shows  that  different  tether 
lengths  can  provide  maximum  protection  if  the  decoy  is  deployed  perpendicular  to  the 
incoming  direction  of  the  missile.  However,  this  is  not  operationally  feasible  under 
practical  flight  conditions.  Nevertheless,  comparison  of  Figure  11  with  Figure  10  shows 
one  thing  in  common.  That  is,  decoy  deployment  has  the  most  affect  on  side-attacking 
missiles  independent  of  the  tether  length.  Actually,  the  longer  the  length  is,  the  better  the 
protection  once  the  reflection  ratio  has  reached  its  effectiveness  equilibrium. 
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Figure  11.  Protection  Envelope  of  Different  Decoy  Deployment  Angles 
for  aJCSj  =  1 
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V.  DECOY  EFFECTS  FOR  AN  EXTENDED  TARGET 


For  most  practical  situations,  the  target  is  not  a  point  reflection  source.  Rather  the 
target  is  extended  and  composed  of  many  scatters.  For  simplification,  the  target  aircraft 
will  be  represented  as  a  10  meters  long  (L  =  10m),  1  meter  radius  (a  =  Im)  cylinder.  Then, 
the  monostatic  reflection  is  given  by  [Ref  1]: 


<T  =  kaL^ 


cosOi 


sin(kLsin  Sj) 
kUsinGi 


2 


(12) 


where  0;  is  the  angle  of  incidence  from  the  broadside  in  a  plane  containing  the  cylinder 
axis,  and  k  =  IkII.,  where  X,  =  0,03m.  This  target  reflection  is  shown  in  Figure  12  where 
CT  =  5.45  m^  and  <&>  =  4.56  m^.  For  the  discussion  in  this  chapter,  the  reflection  ratio 

mtx  tv 

is  defined  as  ai/max(a.j.). 


Figure  12.  Extended  Target  Reflection 
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A.  EFFECTS  OF  REFLECTION  RATIO 


As  shown  in  Figure  13,  the  effects  on  miss  distance  from  a  change  in  reflection 
ratio  achieved  stabilization  when  Op/max(o^)  =  1.5  for  a  side-attack  situation  (missile  from 
90")  and  decoy  at  180".  The  result  is  consistent  with  the  point  target  case.  It  also  shows 
that  when  aD/max(aj)  =  1,  it  protects  the  target  for  each  tether  length.  Thus,  the  reflection 
required  for  the  decoy  in  operation  should  not  be  less  than  that  of  the  target  it  protects. 


B.  CHANGE  OF  MISS  DISTANCE  DUE  TO  DIFFERENT  MISSILE  ATTACK 
ANGLES 

Figure  14  shows  that  for  a  <Tj/max(a^)  ratio  of  1.5  and  the  decoy  at  180",  the 
protection  envelope  provided  by  the  decoy  is  almost  the  same  as  Figure  10.  Still,  it 
provides  maximum  protection  for  a  side-attack  situation,  i.e.,  when  the  decoy  is  deployed 
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perpendicular  to  the  incoming  direction  of  the  missile.  Since  the  result  is  obtained  when 
the  reflection  of  the  decoy  has  reached  its  stabilized  condition,  it  is  safe  to  say  that  the 
towed-decoy  is  a  very  effective  countermeasure  against  a  monopulse  seeker  missile. 


Figure  14.  Protection  Envelope  of  Different  Missile  Attack  Angles 
for  Extended  Target 


C.  CHANGE  OF  MISS  DISTANCE  DUE  TO  DIFFERENT  DECOY 
DEPLOYMENT  ANGLES 

Figure  15  shows  the  protection  envelope  when  <Tj)/max(<Tj)  is  1.5  with  the  missile 
approaching  from  180°.  Compared  with  the  protection  envelope  in  Figure  11,  some  of  the 
notches  have  disappeared  for  longer  tether  lengths.  Also,  the  protection  envelope  grows  as 
the  tether  length  increases.  This  results  in  the  decoy  having  a  wider  deployment  choice 
when  the  tether  length  of  the  decoy  is  longer. 
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Miss  Distance,  meters 


90120 


Solid  -  L=  100m 
Dashed  --  L=  75m 
Dotted  --  L=  50m 


VI.  CONCLUSIONS 


This  thesis  has  investigated  two  kinds  of  target/decoy  reflection  conditions.  One  is 
when  both  the  target  and  decoy  are  point  reflectors,  while  the  other  is  for  an  extended 
target  while  the  decoy  remains  as  a  point  reflector.  Three  different  deployment  cases  are 
considered;  (1)  the  change  of  the  reflection  from  the  decoy,  (2)  the  change  of  direction  of 
release,  and  (3)  the  change  of  initial  missile  incoming  direction,  under  diSerent  tether 
lengths. 

The  result  has  shown  that  a  towed-decoy  is  indeed  an  effective  countermeasure 
against  monopulse  seeker  missiles.  The  operational  parameters  have  the  following 
characteristics; 

•  The  decoy  starts  to  take  effect  when  its  reflection  is  stronger  than  the  protected 
target.  Longer  tether  length  requires  stronger  reflections  to  reach  stabilization 
than  shorter  length.  The  reflection  ratio  required  for  stabilization  is  about  0.1 
for  every  10  meters  length  increment.  When  the  reflection  ratio  reached  1.5, 
the  missile  will  home  directly  onto  the  decoy. 

•  For  both  target  reflection  cases,  the  decoy  is  most  effective  when  it  is  deployed 
toward,  or  perpendicular,  to  the  missile  incoming  direction.  However,  for  the 
extended  target  case,  the  decoy  has  a  wider  choice  of  deployment  regions  than 
the  constant  target  case. 

•  The  decoy  provides  a  very  wide  protection  region  for  any  missile  attack  angle 
except  for  head-on  situations  where  no  protection  is  provided.  There  are 
notches  in  protection  envelopes.  However,  for  longer  tether  lengths,  the 
notches  become  less  pronounced  which  indicates  that  a  longer  tether  length  is 
preferred.  If  we  take  maneuverability  into  account,  tether  lengths  should  be  less 
than  100  meters  but  at  least  twice  as  long  as  the  missile's  lethal  range. 
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Vn.  FUTURE  STUDIES 


Towed-decoys  have  attracted  the  attention  of  many  EW  engineers  and  companies. 
According  to  [Ref  2],  "Towed  decoys  top  the  RF  development  list."  This  thesis  provides 
a  preliminary  insight  regarding  the  operational  deployment  parameters  required  for  the  use 
of  towed  decoys.  There  is  still  additional  work  required  beyond  that  covered  in  this  thesis 
to  ensure  effective  use  of  towed  decoys.  The  following  describes  some  of  this  work: 

•  In  this  thesis  we  assume  a  fiilly  deployed  towed  decoy.  This  results  in  the 
missile  viewing  the  decoy  at  its  fully  stabilized  condition  (i.e.,  fiill  length,  full 
strength  of  reflection).  The  discrimination  conditions  which  result  when  a  single 
source  transitions  to  a  double  source  should  be  further  studied. 

•  The  imaginary  part  induced  (see  Eqs.  5  and  6)  from  the  introduction  of  the 
second  source  can  serve  as  a  warning  that  the  seeker  has  been  jammed.  This 
possible  countermeasure  from  the  seeker's  point  of  view  can  be  studied  as  a 
possible  ECCM  against  towed  decoys. 

•  This  thesis  does  not  consider  the  effects  of  target  glint  and  scintillation.  They 
may  be  major  factors  in  the  determination  of  the  seeker's  resolving  capability 
before  lock-on.  These  effects  need  further  exploration. 

•  Target  maneuvers  which  change  the  relative  bearings  between  the  target  and  the 
decoy  require  further  investigation  to  determine  their  effects  on  the  protection 
envelopes. 
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APPENDIX  A  -  THESIS  MAIN  PROGRAM 


angjnit  =  input('Missile's  Initial  Firing  Angle  = '); 

RO  =  4000; 

RXO  =  5000, 

RYO  =  2000; 

tf  =  input('End  Time  = '); 

md  =  input('End  Game  Requirement  =  ’); 

TL  =  input('To wed-line  Length  = '); 

Dployl  =  input('Decoy  Deployment  Angle  = '); 

RCS  =  inputCRCS  Ratio  (RCSd/RCSt)= '); 

RCSmax=  5  .4562;  %  Target  Maximum  RCS 

RCSd=  RCS  *RCSmax; 

%  Angle  Transformation 
d2r  =  pi/180; 
r2d  =  1 80/pi, 

%  Missile  Initial  Position 

RXinit  =  RX0-R0*cos(ang_init*d2r); 

RYJnit  =  RY0-R0*sin(ang_init*d2r), 

%  Missile  Parameters 
mass  =  1 50; 

Izz  =  1 50; 

PPN  =  4;  %  PPN  Constant 

%  Initial  Run  for  Initial  Conditions 
n=l; 

step(n)  =  .05; 
trun(n)  =  (n-l).*step(n); 
psim(n)=ang_init  *  d2r ; 

VMXI(n)  =  250*cos(psim(n)); 

VMYI(n)  =  250*sin(psim(n)); 

VM(n)  =  sqrt(VMXl(n).^2+VMYI(n).''2); 

RMXI(n)  =  RXjnit; 

RMYl(n)  =  RYJnit; 

RM(n)  =  sqrt(RMXI(n).^2+RMYI(n).^2); 

12B0  =  tm  1  (psim(n)),  %  Transformation  Matrix 
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[VMXB(n),VMYB(n)]  =  i2btran(VMXl(n),VMYl(n),12B0); 
VB(n)  =  sqrt(VMXB(n).^2+VMYB(n).^2), 

AMXB(n)  =  0; 

AMYB(n)  =  0; 

[AMXI(n),AMYI(n)]  =  b2itran(AMXB(n),AMYB(n),I2B0); 
RDOT(n)  =  0, 
r(n)  =  0; 

beta(n)  =  aoa3(VMXB(n),VMYB(n)); 
delr(n)  =  0; 

DR(n)  =  0, 

CFX(n)  =  0; 

CFY(n)  =  0, 

CMZ(n)  =  0; 

FX(n)  =  0, 

FY(n)  =  0; 

N(n)  =  0, 

%  Target  Parameters 
VTXI  =  0, 

VTYI  =  0; 

RTXI(n)  =  RXO, 

RTYI(n)  =  RYO; 

RT(n)  =  sqrt(RTXI(n)/'2+RTYI(n)/^2), 

%  Decoy  Parameters 
dploy  =  dployl .  *d2r, 

TDXI  =  TL.*cos(dploy); 

TDYI  =  TL.*sin(dploy); 

%  Missile/Target  Relative  Geometry 
RMTX(n)  =  RTXI(n)-RMXI(n), 

RMTY(n)  =  RTYI(n)-RMYI(n), 

RMT(n)  =  sqrt(RMTX(n).^2+RMTY(n).^2), 
[mtbx(n),mtby(n)]  =  i2btran(RMTX(n),RMTY(n),I2B0); 
psilos(n)  =  atan2(mtby(n),mtbx(n)); 
psilosd(n)  =  psilos(n)*r2d; 

RCSt(n)  =  raylei(psilos(n)), 

%  Missile/Target/Decoy  Relative  Geometry 
RDXl(n)  =  RTXI(n)+TDXI; 

RDYl(n)  =  RTYl(n)+TDYI, 

MDXI(n)  =  RDXI(n)-RMXI(n), 

MDYI(n)  =  RDYI(n)-RMYl(n); 
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RMD(n)  =  sqrt(MDXI(n).^2+MDYI(n).^2); 

[mdbx(n),mdby(n)]  =  i2btran(MDXI(n),MDYI(n),I2B0); 

psid(n)  =  atan2(mdby(n),mdbx(n)); 

psidd(n)  =  psid(n)*r2d; 

disc(n)  =  max(RMT(n),RMD(n)); 

RXI(n)  =  RMTX(n)+.5.*TDXI; 

RYI(n)  =  RMTY(n)+.5.*TDYI; 

RR(n)  =  sqrt(RXI(n),^2+RYI(n).^2); 

subang(n)  =  acos(  (RMT(n)/^2+RMD(n).^2-TL^2)./(2*RMT(n)*RMD(n)) ); 
dL(n)  =  RMD(n)-RMT(n).*cos(subangl(n)); 
dphase(n)  =  rem(4.  *pi.  *dL(n)./(.03),2*pi); 

VMTX(n)  =  VTXI-VMXI(n); 

VMTY(n)  =  VTYI-VMYI(n); 

%  for  Seeker 
skrang(n)  =  0; 
sd(n)  =  0; 
acmd(n)  =  0; 
ve(n)  =  0; 

Vsum(n)  =  0; 

Vdel(n)  =  0; 

Pest(n)  =  0; 

%  Start  Running 
while  trun(n)  <=  tf 
n=n+l; 
count  =  n; 

step(n)  =  ftime(disc(n-lXTL); 
trun(n)  =  trun(n-l)+step(n); 

[skrang(n),sd(n),ve(n),Vsum(n),Vdel(n),Pest(n)]=yskr(psilos(n-l),skrang(n-l),sd(n-l), ... 

flag,step(n),RCSt(n- 1  ),RCSd,subang(n- 1  ),dphase(n- 1  ),psid(n- 1 )); 
acmd(n)  =  skrang(n)*PPN*VB(n-l), 
delr(n)  =  fin(acmd(n),delr(n-l),count); 

DR(n)  =  delr(n)*r2d; 

[CFX(n),CFY(n),CMZ(n)]  =  aero(beta(n- 1  ),delr(n),  VB(n- 1  ),flag,r(n- 1  ),trun(n)); 
[FX(n),FY(n),N(n)]  =  fin22(CFX(n),CFY(n),CMZ(n),flag,trun(n),  VB(n- 1 )), 

AMXB(n)  =  FX(n)/mass; 

AMYB(n)  =  FY(nymass; 

RDOT(n)  =  N(n)/lzz; 

r(n)  =  intgr(RDOT(n-l),RDOT(n),r(n-l),flag,step(n)); 

I2B  =  tml(psim(n-l)), 

[AMXI(n),AMYI(n)]  =  b2itran(AMXB(n),AMYB(n),I2B), 
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VMXI(n)  =  intgr(AMXI(n-l),AMXI(n),VMXl(n-l),flag,step(n)); 

VMYI(n)  =  intgr(AMYI(n-l),AMYI(n),VMYI(n-l),flag,step(n)); 
psim(n)  =  atan2(VMYl(n),VMXI(n)); 

[VMXB{n),VMYB(n)]  =  i2btran(VMXI(n),VMYI(n),I2B), 

VB(n)  =  sqrt(VMXB(n).^2+VMYB(n).^2), 

RMXI(n)  =  intgr(VMXI(n-l),VMXI(n),RMXI(n-l),flag,step(n)); 

RMYI(n)  =  intgr(VMYI(n4),Wn(n),RMYl(n-lXflag,step(n)), 

RM(n)  =  sqrt(RMXI(n).^2+RMYI(n).^2); 
beta(n)  =  aoa3(VMXB(n),VMYB(n)); 

RTXI(n)  =  RXO+VTXI*trun(n), 

RTYI(n)  =  RYO; 

RT(n)  =  sqrt(RTXI(n)/'2+RTYI(n).^2), 

RMTX(n)  =  RTXI(n)-RMXI(n), 

RMTY(n)  -  RTYI(n)-RMYI(n), 

RMT(n)  =  sqrt(RMTX(n).^2+RMTY(n)/'2); 

[mtbx(n),mtby(n)]  =  i2btran(RMTX(n),RMTY(n),I2B); 
psilos(n)  =  atan2(mtby(n),mtbx(n)), 
psilosd(n)  =  psilos(n)*r2d; 

RCSt(n)  =  rayIei(psilos(n)); 

%  Missile/Target/Decoy  Relative  Geometry 
RDXI(n)  =  RTXI(n)+TDXI; 

RDYI(n)  =  RTYI(n)+TDYl; 

MDXI(n)  =  RDXI(n)-RMXI(n); 

MDYI(n)  =  RDYI(n)-RMYI(n); 

RMD(n)  =  sqrt(MDXI(n).^2+MDYI(n)/^2);  [ 

mdbx(nXmdby(n)]  =  i2btran(MDXl(n),MDYI(n),I2B); 

psid(n)  =  atan2(mdby(n),mdbx(n)), 

psidd(n)  =  psid(n)*r2d; 

disc(n)  =  max(R]^(n),RMD(n)); 

RXI(n)  =  RMTX(n)+.5.*TDXI; 

RYI(n)  =  RMTY(n)+.5.='TDYI; 

RR(n)  =  sqrt(RXI(n)/'2+RYI(n).^2); 

subang(n)  =  acos(  (RMT(n).^2+RMD(n),^2-TL^2)  /(2*RMT(n)*RMD(n)) ); 
dL(n)  =  RMD(n)-RMT(n)  *cos(subang](n)); 
dphase(n)  =  rem(4.*pi.*dL(n)./(.03),2*pi),  %  Wavelength  =  0.03m 

VMTX(n)  =  VTXI-VMXI(n); 

VMTY(n)  =  VTYI-VMYI(n); 

%  End  Game  Condition 

if  ((RMT(n)  <-  md)  |  (RMD(n)  <=  md)) 

MDl  =RMT(n) 

MD2  =  RMD(n) 
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TIME  =  trun(n)  break 


%  end  md  requirement 


end 

end 

end  %  end  while 

min(RMT) 

min(RMD) 
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APPENDIX  B  -  SEEKER  MODEL 


[skrang,skrd,VE,Vsum,Vdel,Pest]==yskr(psit,oldang,rold,flag,step,RCSt,RCSd,... 

subang,dphase,psid) 

r2d=  180/pi; 
d2r  =  pi/180; 
bw  =  .03/.15, 
sq  =  bw/3; 

%  Realtive  to  Seeker's  Tracking  Axis 
psil  =  psit-oldang; 

%  Target  Contributions 
pal  =psil-sq; 

51  =  RCSt.*exp(-2.76,*(pal./bw).^2); 
pa2  =  psil+sq; 

52  =  RCSt.*exp(-2.76.*(pa2./bw).^2); 

%  Realtive  to  Seeker's  Tracking  Axis 
psi2  =  psid-oldang; 

%  Decoy  Contributions 
pbl  =  psi2-sq, 

J1  =  RCSd.*exp(-2.76.*(pbl./bw)/''2).*exp(j.’^dphase); 
pb2  =  psi2+sq; 

J2  =  RCSd.*exp(-2.76.*(pb2./bw).^2).*expO.*dphase), 

%  Antenna  Ouput  Voltages 
VI  =  S1+J1; 

V2  =  S2+J2, 

Vsum  =  V1+V2; 

Vdel  =  Vl-V2; 

%  Error  Voltage 
VEl  =  real(Vdel./Vsum); 

if(abs(VEl))>=  1 

VE  =  0.9987*sign(VEl); 
else 
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end 


VE  =  VE1; 


%  Angle  Estimation 

temp  =  3*bw.*(log(l+VE)-log(l-VE)); 

Pest  =  temp,/(4*2.76); 

%  Seeker  Angle  Servo 
w  =  2*pi*20, 
s  =  exp(-w*step); 
skrd  =  VE*Pest+rold  *s; 
skrangl  =  oldang+Pest; 

%  Angle  Limiter 
if  (abs(skrangl))  >  45*d2r 

skrang  =  45*d2r*sign(  skrangl); 

else 

skrang  =  skrangl; 
end 
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